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Abstract. Within unbroken SUSYQM and for zero factorization energy, I present an 
iterative generalization of Mielnik's isospectral method by employing a Schrodinger true 
zero mode in the first-step general Riccati solution and imposing the physical condition 
of normalization at each iterative step. This procedure leads to a well-defined multiple- 
parameter structure within Mielnik's construction for both zero modes and potentials. 
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The supcrsymmetric procedures arc an interesting and fruitful extension of (one- 
dimensional) quantum mechanics. For recent reviews see pj. These techniques 
are, essentially, factorizations of one-dimensional Schrodinger operators, first dis- 
cussed in the supersymmetric context by Witten in 1981 and well known in the 
mathematical literature in the broader sense of Darboux covariance of Schroedinger 
equations ||. 

In 1984, Mielnik j| introduced a different factorization of the quantum har- 
monic oscillator based on the general Riccati solution. As a result, Mielnik ob- 
tained a one-parameter family of potentials with exactly the same spectrum as that 
of the harmonic oscillator. However, even though in the same year Nieto discussed 
the connection of such a factorization with the inverse scattering approach, and 
Fernandez applied it to the hydrogen atom case, Mielnik's result remained a cu- 
riosity for a decade during which only a few authors paid attention to it. On the 
other hand, constructing families of strictly isospectral potentials is an important 
possibility with many potential applications in physics |l| . This explains the recent 
surge of interest in this supersymmetric issue || . My goal in this work is to give a 
multiple-parameter generalization of Mielnik's procedure based on the ground-state 
function of any soluble one-dimensional quantum mechanical problem. This is just 
a form of Crum's iterations, i.e., repeated Darboux transformations. Some work 
along this line has already been done by Keung et al. Q , who performed an itera- 
tive construction for the reflectionless, solitonic, sech potentials and the attractive 
Coulomb potential presenting relevant plots as well. However, they first go n steps 
away from a given ground state and only afterwards perform the n steps backwards. 
On the other hand, Pappademos et al. 0, working in the continuum part of the 
spectrum, got one- and two-parameter supersymmetric families of potentials strictly 
isospectral with respect to the half-line free particle and Coulomb potentials and 
focused on the supersymmetric bound states in the continuum. Their procedure 
is closer to the method I will present in the following. For more recent works see 
Bagrov and Samsonov || , Fernandez ei aZ [^] , Junker and Roy [^0) , and Rosas-Ortiz 

0- 

In the following, I first briefly recall the mathematical background of Mielnik's 
method and next pass to a simple multiple-parameter generalization for the partic- 
ular but physically relevant zero-mode case. 

I begin with the "fermionic" Riccati (FR) equation y = —y 2 + Vi(x) [the 
"bosonic" one being y = y 2 + Vo(x)] for which I suppose to know a particular 
solution yQ. Notice also that I do not put any free constant in the Riccati equa- 
tions, that is, I work at zero factorization energy. Let us seek the general solution 
in the form y\ = w\ + z/o- By substituting y\ in the FR equation one gets the 
Bernoulli equation —w 1 = w\ + (2yo)w\. Furthermore, using W2 — 1/wi, we obtain 
the simple first-order linear differential equation w 2 — (2yo)w2 — 1 = 0, which can 

be solved by employing the integration factor Fq (x) = e fc 2yo , leading to the so- 
lution W2(x) = (A + f Fo(z)dz) / Fq(x) , where A occurs as an integration constant. 
In applications the lower limit c is either — oo or depending on whether one deals 
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with full- line or half-line problems, respectively. In the latter case, A is restricted 
to be a positive number. Thus, the general FR solution reads 

- J 2 yo p / nx \ 

Vi=ya + r , =y + ° =y + D]n[X+ (F ) , (1) 

where D = It is easy to reach the conclusion that the particular FR solution 
yo corresponds to Witten's superpotential |Q, while the general FR solution y\ = 
yo H ft — is of Mielnik type B. This is especially clear when one is able to 

A+j^ F 

1/2 

identify fo = F with the quantum mechanical ground-state wavcfunction uq of 
the problem at hand. This requires suitable asymptotic behaviour of the Riccati 
solution yo and applying the normalization condition to fo, turning it into a true zero 
mode. As is well known, this case corresponds to the so-called unbroken SUSYQM, 
which will be assumed to hold henceforth. Moreover, — 2y Q (= — lni*^ 2 ) is 
the Darboux transform contribution to the initial Schrodinger potential, i.e., V\ = 
V ~ 2y . Also, the modes 

A + J c F X + J c u 2 

can be normalized and therefore considered as ground-state wavefunctions of the 
bosonic family of potentials corresponding to Mielnik's parametric superpotential. 
The one-parameter true zero modes read 

VAFu /2 ^/Auo f . 

where VA = yf A(A + 1) is the normalization constant. 

Moreover, — 2y l can be thought of as the general Darboux transform contribution 
to the initial potential generating the bosonic strictly isospectral family, which reads 

y« . Va(x) - 2 ^ b (, + /• „ S ) . Va(x) - + . (4) 

This family of potentials can be seen as a continuous deformation of the original 
potential, because the latter is included in the infinite limit of the deforming pa- 
rameter A and v± oa = uq as well. In more intuitive terms, Mielnik's method based 
on an intial Schroedinger true zero mode may be called a double Darboux technique 
of deleting followed by reinstating a nodeless ground-state wavefunction uq(x) of 
a potential Vo(x) by means of which one can generate a one-parameter family of 
isospectral potentials V\(x), where A is a labeling, real parameter of each member 
potential in the set. 

One can go on with one of the strictly isospectral bosonic zero modes u\ 1 = 
"/i, (i.e., by choosing A = Ai) and repeat the strictly isospectral construction, 

Al+ Je "O 
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getting a new two-parameter zero mode 

u Ai,A 2 = 



U 



A 2 + /> 2 Al (A 1 + />g)(A 2 + /;<) 
The two-parameter true zero modes read 



VAlA 2 U 



(Ax+r^KAa+r^) • 

The resulting two-parameter family of strictly isospectral potentials will be 



V Xu x 2 =V -2—\n 



Ai 



A 2 



(5) 



(6) 



(7) 



At the ith -parameter level, one will have 



^ A a = Vo-2^1n 



and 



Ai+ (a 2 + f<)- (' 
VAiA 2 ...A t M 



A* + / ^Ai...A«_i 



(8) 
(9) 



Explicit formulas for the parametric zero modes can be obtained if one uses a 
notation based on the integration factor f a F = T{x) — T{c) = A X T . Then 



v\A x ) = 



/ A7?i 



Ai + A X T ' 



(10) 



Next, one can calculate 



(A-^(c)+^(x)) 2 



1 



A X T 



Thus 



WAlA^x) 



1 

A ~ A + A X T A(A + A x jT) 
VAiA 2 u 



A1A2 + (Ai + A 2 + 1)A X JF 



At the next step one gets 



VAiA 2 A 3 u 



AiA 2 A 3 + (AiA 2 + A 2 A 3 + A 3 Ai + Ai + A 2 + A 3 + l)A x T 
and the general formula at the i level can be written down in the form 



V\ u \2,...,\i(x) 



\ZAi....AjM 

C^+C^A^ 



(11) 
(12) 

(13) 
(14) 



4 



where the first coefficient in the denominator is the product of all parameters, 
whereas the second coefficient is just the sum over all the rest of lower order Viete- 
type expressions in the parameters. By the same token, one can write a general 
formula for the strictly isospectral potentials 



\\..v, a = Vo-2D 2 ln[C^+C ( 2 l >A x f] = V Q 



Cf ] + cfA x T (C^ + ci l) A x f) 2 



= v Q - ,^ u ^ + M£ 



(15) 

which may be considered as the generalization of the furthest right-hand side of 
Eq. (4) and represents a simple generalization of Mielnik's one-parameter potentials. 
Since Ai...Aj = C^(C^ + Cf), one might think that there is nothing new in 

(14) and (15) with respect to a common Mielnik solution with an effective parameter 
= jcf. However, I will argue that by performing such an equivalence one 

loses a certain type of information. This information is a consequence of the sym- 
metry of (14) and (15) in the space of parameters. One can see that the subindiccs 
of any pair of parameters can be interchanged without affecting the formulas. Thus, 
each A parameter can be varied independently of the others, making it possible to 
put questions related to the following type of situation. Suppose we construct two 
Mielnik potentials corresponding to Ai and A2 and ask what is the potential bearing 
true zero modes that for Ai — > ±00 goes to the Mielnik case for A2, whereas for 
A 2 — > ±00 it goes to the Mielnik case for Ai . The answer is provided by the construc- 
tion of this work and corresponds to the particular case V\ lt \ 2 bearing the true zero 
modes w\i,A 2 - Indeed, as one can easily check, VAi,±oo — ^Ai and V\ x .±oa = V\f > 
whereas w±oo,a 2 = v x 2 an d V±oo,\ 2 = V££ . In the general case, one starts with a set 
of i Mielnik potentials corresponding to i fixed values of Mielnik's parameter and 
asks the same question, this time for the set of i asymptotic limits. The answer is 
given by (14) and (15) and cannot be provided if one works with only one effective 
parameter unless its multiple-parameter value found above is used. 

Another interesting remark is the one-to-one relationship between any polyno- 
mial equation a^x 1 + aix' 1 ^ 1 + ... + a» = and the present iterative construction. 
If we consider the A parameters as the zeros of such arbitrary polynomials, we can 
write (14) as 

(-l)X + (}J (-l) l ai)A x .F 

and in (15) one can substitute the same type of denominator as in (16). There is 
only one constraint on the employed polynomials, which one should impose in order 
to avoid possible singularities. Usually the integral A X T in the denominators of 

(15) and (16) is of the kink type, i.e., it may be written in the form a + (3K(x), 
where the function K(x) has the kink behavior, taking values between -1 and +1 
and a and j3 are some constants, of which a may be zero. Then the allowed intervals 
for the effective parameter arc X^ ff > /3 - a and \ f ff <-((3 + a). When a = 

one gets |A^| > (3. 
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It is worthwhile to mention that the previous iteration process can be understood 
most easily from the Riccati equation standpoint as follows. To get, for example, 
the two-parameter zero mode, one should start again with the FR equation y = 
—y 2 + V\{x) and take as the known particular solution y^ = yo + y\ t , where 
y\, = 78 • The intermediate Bernoulli equation will be —w 1 = w\ + 2y^wi. 

Xi+ J (Fo) 

This is turned into a first-order differential equation by the inverse function method. 

(_ f x 2-»< 1 ») 

The integration factor of the latter is = cxp J <= p and the solution for the 
first-order differential equation can be written u>2 = (A2 + f x F\ 1 dz)/F\ 1 . From 
this presentation it is clear how one should proceed for an arbitrary step. Also, 
the logarithmic derivative notation in Eq. (1) is equally convenient to have a clear 
image of the iteration process. Thus, one can generate hierarchies of parametric 
Schrodinger zero modes of any desired order by means of the Riccati connection. 

The parametric normalization deletes the interval [—1,0] from the parameter 
space of ^gjf At the —1 limit, one can make a connection with the Abraham- 
Moses isospectral technique 113] , whereas at the limit the connection can be done 
with another isospectral construction developed by Pursey []l3| . This connection is 
only from the point of view of the potentials; the zero modes as worked out here 
just disappear. 

In conclusion, I have shown explicitly the way Crum's iteration works when the 
general Riccati solutions (general superpotentials) at zero factorization energy are 
based on the corresponding Schrodinger ground-state wavefunctions, obtaining gen- 
eral formulas for this simple 'generalization' of Miclnik's one-parameter SUSYQM 
isospectrality. Plots of the two-parameter formulas for the harmonic oscillator case 
are presented in Fig. 1-5. One may consider the results of this work as pointing to 
an interesting hierarchical structure within the general Riccati solution produced 
by a particular type of repeated Darboux transformations when the normalization 
condition of quantum mechanics is taken care of at each iterative step. 
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Fig. 1 

Two-parameter strictly isospectral harmonic oscillator potentials (h = m = u> = 1), for 
fixed A2 = 0.2 and Ai € [0.1 , 5]. They are identical to Mielnik harmonic oscillator 
potentials with X ( ^ f e [0.0154,0.1613]. 
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Fig. 2 

The corresponding true zero modes. 
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Fig. 3 

The two-parameter true zero modes at fixed x = —1.4 as a function of the two 

parameters. 
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Fig. 4 

Same modes as in Fig. 3 at fixed x = —1.6. 
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